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LIGHT  SCATTERING  FROM  A  DEEP  METALLIC  GRATING 
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Thomas  F.  George 

Departments  of  Chemistry  and  Physics,  State  University  of  New  York 
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ABSTRACT 

The  conditions  under  which  the  Rayleigh  hypothesis  is  exact 
and  the  convergence  properties  of  the  Rayleigh  expansion  are 
considered.  The  identification  of  the  cause  of  the  deficiency  of 
this  expansion  suggests  an  alternative  "dressed"  expansion  with 
presumably  simpler  convergence  properties.  This  proposition  is 
checked  for  a  sinusoidal  grating  (SG),  for  which  convergence  is 
found  for  an  arbitrary  value  of  B  *  2iTg/d,  where  g  and  d  denote  the 
height  and  periodicity  of  the  SG,  respectively.  The  dressed 
expansion  is  used  to  analyze  the  surface  plasmon  dispersion  and 
local  field  enhancement  distribution  pertaining  to  the  SG  in  the 
limit  as  B  goes  to  infinity.  The  dispersion  relation  is  comprised 
of  two  bands.  The  local  field  enhancement  predicts  stronger  fields 
at  the  bottoms  of  the  troughs  than  at  the  peaks  of  the  SG. 

INTRODUCTION 

Light  scattering  from  a  grating  is  a  well -developed  subject  in 
classical  optics.  The  underlining,  well-known  physical  process  is 
(elastic)  Bragg  scattering,  i.e.,  the  exchange  of  surface-parallel 
momentum  between  the  grating  and  the  incident  light  by  an  amount 
nk_,  where  n  *  0,il,i2,...  and  the  grating  momentum  is  k^  «  2ir/d. 
All  notations  are  defined  in  Fig.  1.  It  is,  however,  important  to 


0 


N 


Fig.  1.  Generic  grating 
defining  the  notations.  The 
incident  light  is  directed 
downward,  and  the  hatched  area 
indicates  the  metal.  The 
domains  "a"-"d"  pertain  to  the 
discussion  of  the  Rayleigh 
expansion. 
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realize  that,  barring  a  few  recent  studies,^  most  of  the  work 
pertains  to  the  restricted  physical  domain  of  shallow  gratings 
where  B  *  2iTg/d  <  1.  Light  scattering  from  such  gratings  is 


qualitatively  different  from  that  pertaining  to  deep  gratings, 
i.e.,  when  3  »  !•  In  the  former  case,  the  grating  can  exchange 
only  a  limited  number  of  quanta,  i.e.,  n  -0(1).  Consequently, 
the  scattered  light  is  comprised  of  a  dominant  specular  and  a  few 
Bragg  components.  On  the  other  hand,  for  deep  gratings,  n  -  9(3) 
»  1  (see  below).  Hence,  the  scattered  light  is  comprised  of  many, 
equally  important  and  strongly  interfering  Bragg  components,  which 
gives  the  scattered  light  a  new  qualitative  character.  The 
theoretical  description  of  the  deep  grating  domain  is  at  the  focus 
of  this  paper. 

RAYLEIGH  EXPANSION  AND  HYPOTHESIS 


The  Rayleigh  expansion^'  provides  a  framework  to  the  shallow 
grating  work  and  a  starting  point  to  our  analysis.  It  results  from 
the  underlying  symmetry  of  the  grating:  Since  it  is  invariant 
under  translations  x  ■*  x-Hid,  where  n  is  any  integer,  the  Floquet- 
Bloch  theorem  implies  that  the  electromagnetic  field  ?(x,z)  has  the 
property  ?(x-fd,z)  >  exp[ik||^]^(x,z)  where  kn  is  a  surface-parallel 
(Fig.  1)  momentum  label.  Consequently,  the  Fourier  series 
expansion  of  the  fields  in  a  "unit  cell"  0  <  x  <  d  (the  z- 

dependence  is  determined  from  the  wave  equation)  suffices  to 
describe  the  field  throughout  the  xz-plane.  When  the  proper 
boundary  conditions  at  infinity  and  the  vector  character  are 
incorporated,  the  ensuing  p-wave  expansion  is  the  Rayleigh 
expansion: 

^  r  i(k,x  -  W  (l)z) 

2  (x,z)  -  )  {C  (1)5  (t)e  ^  “ 

t»-«» 

i(k  X  +  W  (t)zl 

+  A  (t)p„  .(l)e  ‘  ®  )  , 

®  o»+ 


where 


k(a)  ■  ’^Ejjk  ,  k 

W  (O  -  [k^(a)-k?]^'2 
a  l 


3  -  gkG, 


W^(t)  -  [k‘(a)-k,]*'‘  ,  k  -  k  +  Ik-  ,  Re(W  ),Im(W  )  >  0, 

(kj^z  ?  W^(t)x],  s  -  X  X  z  .  (2) 

In  (2),  X  and  z  are  unit  vectors  in  the  x-  and  z-  directions 
(Fig.l),  and  a  denotes  a  domain  in  the  xz-plane  with  a  constant 
dielectric  constant  c  .  Once  a  convenient  partition  of  the  xz- 
plane  into  domains  has  been  chosen,  the  coefficients  C  ,A  are 
determined  by  matching  the  boundary  conditions  across  the  domain 
boundaries.  To  avoid  mathematical  pitfalls  in  the  subsequent 
discussion,  we  confine  ourselves  henceforth  to  gratings  such  that 
the  fields  are  non-singular  everywhere. 

The  Rayleigh  hypothesis  can  now  be  introduced  in  terms  of 
expansions  (1).  For  example,  the  grating  in  Fig.  1  calls  for  four 
expansions  (1)  while  the  Rayleigh  hypothesis  asserts  that  only  two 


t  ^ 


are  needed  since  the  expansions  in  domains  "a"  and  "b"  are 
identical  and  likewise  for  domains  "c”  and  "d".  By  matching  the 
boundary  conditions  across  the  (z*g)-plane,  we  conclude  that  the 
hypothesis  is  exact  for  gratings  whose  profile  is  expressible  in 
terms  of  a  finite  Fourier  series.  For  other  gratings,  the  Rayleigh 
hypothesis  may,  or  may  not  be  exact. 

DRESSED  RAYLEIGH  EXPANSION 

Attempts  to  apply  the  Rayleigh  expansion  to  deep-grating 
calculations  inidicate  its  divergence  for  3  >0.6  (sinusoidal 
gratings).  This  si^u^j^ion  has  motivated  the  introduction  of 
alternative  schemes^  *  in  the  quest  to  enlarge  the  R-range 
accessible  for  calculations.  Our  approach  is  to  expose  the  origin 
of  the  instability  of  (1)  when  3  »  1,  which  then  naturally  leads 
to  an  alternative  expansion  the  dressed  expansion  --  with 
excellent  convergence  properties. 

For  the  grating  in  Fig.  1,  the  expansion  in  domain  "0",  which 
combines  domains  "a"  and  "b",  is 


alternative  expansion  -Jq  the  dressed  expansion  --  with 


2g(x,z) 


i(k^x  + 


Wp(t)z] 


Disregarding  the  uninteresting  phase  exp(ik,x)  and  geometrical 
factor  Pq  ^(t)  in  (3),  the  |t|  -►  •  behavior  is  Csee  (2)): 

iW.(t)z  -k_|t|z 

A_(t)e  ^  A-(l)e  .  (4) 

”  |ti-  ” 


iWQ(t)z 


Consequently,  at  the  bottom  of  a  trough,  where  ~g  <  z  <  0,  the  z- 
dependent  factor  in  (4)  diverges  exponentially.  Therefore,  to 
render  the  total  field  Eq(x,z)  finite,  the  exact  A^Ct)  must 
converge  at  least  as  exp(-3|t|).  Hence  for  z  <  0,(3)  is  a  sum  of 
energy  terms,  most  of  which  are  products  of  exponentially  large 
factors  times  exponentially  small  numbers.  The  calculated  A.(£), 
however,  always  entail  an  error  which,  when  multiplied  by 
exp(3|l>|)t  yields  large  erros  in  EQ(x,e).  Thus,  (3)  is 

I  scheme  for 


intrinsically  unstable  as  a  scheme  for  3  »  1  (man 

calculations. 

This  deficiency  is  easily  remedied  by  rewriting  (4)  as 


(many  |i(' 


AQ(l)e 


1Wq(1)z 


aQ(t)e 


lWQ(l)(z+g) 


where  oiqCI’)  “  exp(“iW  (t)z}.  By  construction  the  RHS  of  (5) 
always  converges  for  z-fg  >  b,  and  when  z-t-g  ■  0,  the  QQ(f)  converge 
to  render,  by  assumption,  a  finite  ]Sq(x,z)  on  the  grating's 
surface.  Therefore,  by  transcribing  vl)  into  the  "dressed" 
expansion,  as  defined  by  the  transformation  (5),  a  major  flaw  in 
(1)  is  eliminated. 


construction  the  RHS  of  (5) 
z+g  ■  0,  the  aj,(£.)  converge 


*  A  Va  -N  -  1  .  . 


The  convergence  of  the  dressed  expansion  ^gn  be  explicitly 
demonstrated  for  a  sinusoidal  metallic  grating.  Lack  of  space 
does  not  allow  us  to  outline  the  analysis.  The  results  are  that 
the  dressed  expansion  converges  for  an  arbitrary  large  3,  and  the 
number  of  significant  components  in  the  dressed  expansion  is  on  the 
order  N  -  43.  This  analysis  explains  why  the  (bare)  Rayleigh 
expansion  diverges  for  3  ^  0.6. 

SURFACE  PLASMON  (SP)  DISPERSION  RELATION  IN  THE  3  "  LIMIT 

As  an  applicati^^  of  the  dressed  expansion,  we  consider  the  SP 
dispersion  relation^  in  the  3  limit  for  ^^sinusoidal  grating. 
Lack  of  space  allows  us  only  quote  the  result: 


3[k^(0)  +  k^(l)}^ 


sin(L(a+y) Jsin[^^]  -  sin(L(a*Y)]slnt^^I  , 

“  { - ^ ,  (6) 

sln(L(o+Y)]sin(^^]  +  sin(L(o“Y)]sin(-^] 

where  k(i)  ■  k/e^^  from  Eq.  (1),  y  ■  2vk/kg  , 

a  -  -~{-(k^(0)  +  k^(l)]}^^^  ,  (7) 

and  L  -  Q(3)  »  1»  the  precise  value  of  which  is  immaterial. 

The  analysis  of  (6)  is  straightforward:  The  LHS  is  a  smooth 
function  of  x  >  k/k  ,  where  k  *  w  /c  and  u>  is  the  (volume) 
plasmon  frequency.  At^x  «  0  it  is^unity?  decreasing  monotonically 
to  -<»  at  X  *_l//2  and  subsequently  increasing  and  becoming  positive 
for  X  >  /3/2.  The  RHS  of  (6),  on  the  other  hand,  is  rapidly 

oscillating  with  frequency  -1/3  and  amplitude  -3>  Furthermore, 
since  a  (Eq.  (7))  is  real  for  x  <  1//2  and  imaginary  for  x  >  I//2, 
the  RHS  of  (6)  is  positive  and  negative,  respectively. 
Consequently,  (6)  has  a_  large  number  of  solutions  (&(3))  in  the 
domains  0  <  x  _<  1/2  and  /2/2  <  x  <  /3/2.  _  The  frequency  band  gap 

between  these  two  band  is  hence  Au/u  ■  (/2~l^j^2.  This  outcome  is 
in  keeping  with  the  known  behavior  fo?  3  «  !• 

In  summary,  we  have  exposed  the  reason  for  the  poor 
convergence  of  the  Rayleigh  expansion  when  3  »  1 i  and  introduced 
an  alternative  expansion  --  the  dressed  expansion  --  which  has 
allegedly  very  good  convergence  properties.  This  premise  is 
explicitly  demonstrated  for  a  sinusoidal  grating.  We  have  also 
discussed  the  surface  plasmon  dispersion  relation  in  the  3  -»  «> 
limit,  where  we  fdnd  that_the  branches  tend  to  cluster  into  two 
bands  separated  by  Au/u  ■  (i/2*‘l)/2. 
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